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Abstract. We report a ~38 ks X-ray observation of McNeil's Nebula obtained with XMM-Newton on 2004 April 4. 
V1647 Ori, the young star in outburst illuminating McNeil's Nebula, is detected with XMM-Newton and appears 
variable in X-rays. We investigate the hardness ratio variability and time variations of the event energy distribution 
with quantile analysis, and show that the large increase of the count rate from V1647 0ri observed during the 
second half of the observation is not associated with any large plasma temperature variations as for typical X- 
ray flares from young low-mass stars. X-ray spectral fitting shows that the bulk (~ 75%) of the intrinsic X-ray 
emission in the 0.5-8 keV energy band comes from a soft plasma component, with fcT so ft = 0.9 keV (0.7-1.1 keV, 
at the 90% confidence limit), reminiscent of the X-ray spectrum of the classical T Tauri star TW Hya, for which 
X-ray emission is believed to be generated by an accretion shock onto the photosphere of a low-mass star. The 
hard plasma component, with fcTh ar d = 4.2 keV (3.0-6.5 keV), contributes ~ 25% of the total X-ray emission, 
and can be understood only in the framework of plasma heating sustained by magnetic reconnection events. We 
find a hydrogen column density of Nn = 4.1 x 10 22 cm -2 (3.5-4.7 x 10 22 cm -2 ), which points out a significant 
excess of hydrogen column density compared to the value derived from optical/IR observations, consistent with 
the picture of the rise of a wind/jet unveiled from ground optical spectroscopy. The X-ray flux observed with 
XMM-Newton ranges from roughly the flux observed by Chandra on 2004 March 22 (i.e. ~10 times greater than 
the pre-outburst X-ray flux) to a value two times greater than that caught by Chandra on 2004 March 7 (i.e. 
~200 times greater than the pre-outburst X-ray flux). The X-ray variability of V1647 0ri in outburst is clearly 
enhanced. We have investigated the possibility that V1647 0ri displays a periodic variation in X-ray brightness as 
suggested by the combined Chandra+ XMM-Newton data set. Assuming that the X-ray flux density is periodic, 
the folding of the two Chandra observed X-ray flux densities with the XMM-Newton ones leads to three periodic 
X-ray light curve solutions. Our best period candidate is 0.72 day, which corresponds to the time scale of the 
Keplerian rotation at a distance of 1 and 1.4 stellar radius for a one solar mass star aged of 0.5 and lMyrs, 
respectively. We propose that the emission measure, i.e. the observed X-ray flux, is modulated by the Keplerian 
rotation of the inner part of the V1647 0ri accretion disk. 

Key words. Stars: individual: V1647 0ri - Stars: pre-main sequence - X-rays: stars 



1. Introduction 

At the end of January 2004, a new 
nebula was discovered serendipitously 
tween M 78/NGC 2068 and the reg 
These two well known star-formin| 
L1630 dark cloud are located in the 
the Orion B giant molecular cloud, 



bright fan-shaped 
<|McNeil 2004|) be- 
ion of HH 24-26. 
regions of the 
northern part of 
at a distance of 
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~400pc ( |Anthony-Twarog 1982] ). At the apex of McNeil's 
Nebula lies the young stellar objec t (YSO) IRA S 05436- 
0007=2MASSJ 05461313-000604 8 llClark 1991|l, associ- 
ated with a faint /-band source ( Eisloff eT fc Mundt"l 997) 
detected in the (sub)millimeter (1.3 mm source LMZ12 
of ILis et al. 19 99, and 850 fim source OriBsmm55 of 
Mitchell et al. 20010 . This YSO is now displaying a dra- 
matic optical/IR outburst, which is the origin of the rise 
of this reflection nebula. Briceho et al. (2004) have con- 
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strained the start of this outburst to the beginning of 
November 2003, and have obtained a light curve in the 
/-band of V1647 0ri (designation of the variable star 
illuminating McNeil's Nebula; ISamus 2004fl showing a 
~5mag brightening in about 4 months. They show that 
the timescale for the nebula to develop is consistent with 
the light-travel time, indicating that we are observing 
light from the central source scattered by material in the 
cometary nebula. 

The nature of the V1647 0ri outburst and its con- 
nection with other pre-main sequence (PMS) eruptive 
objects, namely EXors (stellar prototype: EXLupi; 
|McLaughlin 1946||Herbig et al. 2001), and FUor s (stellar 
prototype: FU Ori~ |Herbig 1966| |Herbig 1977), is still 
debated 
and 



at the moment. 



McGehee et al. (2004) 



Reipurth fc Aspin (2004) 



have noted a resem- 



blance to EXors: whereas 



Abraham et al. (2004) 



Kun et al. (2004) have 



Briceno et al. (2004) 
Walter et al. (2004)[ and 
proposed a FUor event. 
Vacca et al. (2004) have even reported that the near-IR 



spectrum of V1647 0ri does not appear similar to any 



known FUor or EXor object. Aspin & Reipurth (2004) 



reported an i?-band observation of V1647 0ri with the 
Gemini telescope on Mauna Kea. These latest data, 
obtained on 2004 August 29 (shortly after the object 
emerged from conjunction with the Sun), indicate that 
V1647 Ori continues to be in a state of elevated optical/IR 
emission, and that McNeil's nebula remains bright, and 
therefore that the duration of the outburst so far exceeds 
about nine months. Indeed, continued optical/near-IR 
monitoring arc necessary to determine the duration of 
the outburst, because outbursts last only several months 
in EXors compared to several decades in FUors. Both 
types of outburst are thought to be driven by a sudden 
increase of accretion through a circumstellar disk (e.g., 
|Hartmann fe Kenyon 19961 an d references therein), but 
the distinction between FUors and EXors is still entirely 
empirical. The interpretation that the outburst event 
has its origin in accretion processes is supported by 
high-resolution IR spectra of V1647 0ri showing CO 
emission lines, likely originating from ~ 2500 K gas in an 
inner accretion disk region where substantial clearing of 
dust has occurred ( |Rettig et al. 2 005). 

The pre-outburst spectral energy distribution (SED) 
of V164 7Qri from IR to milli meter shows a flat-spectrum 
source \ Abraham et al. 2004 1. This kind of flat SED is 
usually interpreted in terms of a circumstellar envelope 
(Ke nyon fc Hartmann 199T] ). Andrews et al. (2004) pro- 



pose that V1647 Ori is a transition object between a proto- 
star with circumstellar disk plus remnant circumstellar en- 
velope (Class I protostar), and a classical T Tauri star with 
circumstellar disk (Class II sources). Spitzer observations 
of V1647 0ri in early March 2004 HMuzerolle et al. 2005|> 
show a factor of 15-20 increase in brightness across the 
spectrum from the optical to 70 /im, leading to a bolo- 
metric luminosity of 44 Lq, i.e. ~15 times higher than the 
pre-outburst level. 



X-ray imaging spectroscopy observations with 
Chandra have revealed a factor ~50 increase in the 
X-ray count rate from V1647 0ri during its outburst 
compared to the pre-outburst state ((Kastner et al. 20 04) . 
The coincidence of a surge in X-ray brightness with 
the optical/IR outburst demonstrates that strongly 
enhanced high energy emission from V1647 Ori occurs 
as a consequence of high accretion rates. The burst of 
X-rays was most probably generated via star-disk mag- 
netic reconnection events that occurred in conjunction 
with such mass infall. This process may also launch 
new, collimated outflows or jets HGoodson et al. 19 97). 
Indeed, before its recent eruption, V1647 0ri had been 
identified as the exciting source of a chain of extended 
emission nebulosity that appears to terminate at HH23, 
a shock-excited Herbig-Haro object, located ~ 3' North 
from V1647Qri ((Eisloftel fc Mundt 19971 ILis et al. 19991 
|Reipurth fc Aspin 2004| ) . The presence of these structures 
suggests that the present optical/IR/X-ray outburst of 
V1647 0ri may be merely the latest of a series of such 
events. Another outburst may have occurred about 
37yrs before the present event, based on a photograph 
obtained in Oct. 1966 which shows a similar cometary 
nebula (Mallas & Krcime r 1978(1 . No clear evidence for 
the presence of a molecular outflow has been found in 
the submillimeter CO spectral line maps l|Lis et al. 1 999 
lAndrews et al. 2004(1 . However the Ha line, which has 
been detected in strong emission, displays a pronounced 
P Cygni profile, with an absorption trough reaching 
velocities up to 600kms _1 ( |Reipurth fc Aspin 2004| ), 
implying significant mass loss in a powerful wind. 

We report here an X-ray observation of V16470ri in 
outburst obtained with XMM-Newton, showing enhanced 
X-ray variability from V1647 0ri. We present in j0 the 
XMM-Newton observation and data reduction, we study 
the variability in and the X-ray spectrum in 21 We 
compare the X-ray flux observed with Chandra and XMM- 
Newton in §5\ and we discuss in ^the origin of the X-ray 
flux variations observed from V1647 Ori. 

2. XMM-Newton observations and data reduction 

We observed V16 47 Ori with XMM- 
Newton l|Jansen et al. 2001(1 on 2004 April 4. The point- 
ing nominal J2000 coordinates were aj2ooo=05 h 46 m 13!l, 
< 5.T2ooo=-00°06'04'.'6. We used the full frame science mode 
of the EPIC cameras with medium optical blocking filter. 
The total exposure times were ~37.0ks and ~38.7ks 
for the EPIC/PN (|Striider et al. 2001(1 and the two 
EPIC/MOS spectro- imagers ((Turner et al. 200~T)l . respec- 
tively. We also acquired optical and UV data with the 
Optical Monitor (OM) with the following filter sequence: 
V, U, B, and UVW1; however the OM data are severely 
contaminated by stray light ((Mason et al. 200~T)l from 
the extended reflection nebula M78, which is located 
~10'-~15' North-East of V16470ri, preventing us from 
using these data at the position of V16470ri. We will 
focus here only on the X-ray detection of V1647 Ori. 
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The data were reduced using the XMM-Newton 
Science Analysing System (SAS, version 6.0.0) 1 . The event 
lists were produced using tasks epchain and emchain. 
The observation is mainly affected by a background pro- 
ton flare having a duration of one hour, centered at 
lh U.T. on 2004 April 4, and peaking on PN 2 to 50- 
150 count ks^ 1 arcmin -2 . Removal of time intervals with 
flaring background above 10 and 5 count ks _1 arcmin~ 2 
for PN and MOS, respectively, leads to effective expo- 
sure times of ~33ks and ~38ks, respectively. X-ray im- 
ages were then produced by selecting good events within 
the fields-of-view of the detectors 3 . 

Vignetting corrections were made from exposure 
maps computed on the corresponding energy bands. 
Figure ^ shows the tricolour X-ray image of the 
region surrounding McNeil's Nebula (see for compar- 
ison Fig. 2 of IKastner et al. 2004). The X-ray source 
XMM J054613. 1-000604 found by the pipeline is lo- 
cated at a J2 ooo=05 h 46 m 13?2, <5j 20 oo=-00 06'04"l with 
a la error box of 0.2", and displays a count rate of 
0.036 ±0.001 counts^ 1 in PN and 0.012 ± 0.001 count s" 1 
in each MOS detector within the energy band 0.2-12 keV. 
This X-ray source is located at 0.9" from the 2MASS 
position of V1647 0ri. Computing the root-mean-square 
dispersion of distances between the positions of X- 
ray sources found by the pipeline and their 2MASS 
counterpart positions, we found the same offset value, 
therefore we will consider this XMM-Newton source to 
be the same X-ray counterpart of V1647 0ri already 
detected by Chandra (Kast ner et al. 20 04) . The other 
X-ray sources visible in Fig. ^ ar e (from left to right): 
XMMJ054619.4-000520, XMM J054618. 9-000539, and 
XMM J054611. 6-000627, which are the X-ray coun- 
terparts to SSV64=LkHa301, GSC2 S020001258, 
and 2MASS J054 6116 2-0006279, respe ctively (see 
IKastner et al. 20041 and lSimon et al. 2004jl . 

We selected in each camera X-ray events of the 
V1647 Ori X-ray source for variability and spectral analy- 
sis. For PN, we selected only single and double pixel events 
(i.e. PATTERN in the to 4 range) with FLAG value equal to 
zero. Source plus background events were extracted within 
a circular region using the SAS task eregionanalyse. The 
extraction position and radius (r = 23") were optimized 
by this task to maximize the signal-to-noise ratio, taking 
for the estimation of the background an annular region 
(?' = 23"-44") centered on this source and located on 
the same CCD, where areas illuminated by other X-ray 
sources were excluded. Background events of this annular 
region were also extracted. 



For a description of the standard procedures see, e.g., the 
XMM-Newton ABC Guide, Snowden et al. (2004), available at 
http : //heasarc . gsf c .nasa.gov/docs/xmm/abc/ . 

1 Background value obtained from the detector background 
light curve provided for PN by epchain. 

3 We selected single, double, triple, and quadruple pixel 
events (PATTERN in the to 12 range) for both PN and MOS, 
plus the FLAG list #XMMEA_EM for MOS. 




05 B 46 m 20' 18' 16" 14" 12" 

"J2000 



Fig. 1. XMM-JVewton/EPIC image of the region sur- 
rounding McNeil's Nebula. The three images obtained 
with PN, MOS1, and MOS2 EPIC cameras on 2004 April 
4 were added after correction from vignetting. The size 
of the image pixel is 3". This image has red, green, and 
blue colour coding for X-ray photons in the 0.5-1.5, 1.5- 
2.4, and 2.4-8.0 keV energy bands, respectively (see algo- 
rithm of Lupton et al. 2004). For comparison purposes 
we overlay a i?-band contour map of the nebula obtained 
with VLT/FORS2 on 2004 February 18. X-ray sources and 
2MASS sources are marked with 'x' and '+' respectively. 
A bright embedded X-ray source is spatially coincident 
with V1647 Ori at the apex of McNeil's nebula. 



3. Variability and quantile analysis 

We limited our variability analysis to the PN observ- 
ing time interval, where all three EPIC cameras were 
working together, which allowed us to build directly an 
EPIC (PN+MOS1+MOS2) light curve. Using the SAS 
task evselect, we first built two light curves with 1 s time 
bins from the source+background and the background 
event lists (where we kept flaring background time inter- 
vals) . These light curves begin exactly at the start of the 
first Good Time Interval (GTI) of the PN aimpoint CCD 
#4 (where the source is detected). We rebinned the ex- 
tracted light curve to a longer time interval with the ftools 
task lcurve to increase the signal. Then, we subtracted 
from the source+background light curve the background 
light curve scaled to the same source extraction area, us- 
ing the ftools task lcmath. We also corrected this source 
light curve for the observing time lost, for example, due to 
the triggering of counting mode during high flaring back- 
ground periods, when the count rate exceeded the detec- 
tor telemetry limit. From the GTI extension, we computed 
with an IDL routine for each time bin the ratio between 
the time bin length and the observing time lost in this 
time bin, and we multiplied count rate and error by this 
linear correction factor. Finally, light curves of the three 
detectors were summed to produce the EPIC light curve. 

The upper panel of Figure [2] shows the EPIC X-ray 
light curve of V1647 Ori in the energy band from 0.5 to 
8.0 keV binned to 1000 s. The source is variable, with a 
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Hour since April 4, 2004 Oh U.T. 

Fig. 2. XMM-Newton/EPIG background subtracted X- 
ray light curves of V1647 Ori. The upper panel shows 
the EPIC (PN+MOS1+MOS2) X-ray light curve of 
V1647 Ori with one sigma error bars in the energy band 
from 0.5 to 8.0 keV. The grey stripes indicate the small 
observational time intervals before the start and the end 
of the PN exposure where only the two MOS detectors 
were observing. The middle panels show the soft (5=0.5- 
2.8 keV) and hard (iJ=2.8-8keV) band X-ray light curves. 
The lower panel shows the variation of the corresponding 
hardness ratio. The dashed line indicates the hardness ra- 
tio for the second half of the observation. 



low level in count rate during the first half of the observa- 
tion, varying between ^0.025 and ~0. 005 counts -1 (note 
that the decrease of count rate around 1 h U.T. cannot 
be a consequence of the high flaring background because 
we have corrected for the losses of observing time due to 
triggering of counting mode). This low level was then fol- 
lowed by a much higher level of count rate during the 
second half of the observation. This latter phase consisted 
of a slow (~4h) rise with an abrupt jump to a peak at 
~0. 13 counts -1 , and then a faster (~lh) decay from the 
peak of ~0.13 count s -1 back down to the general level of 
enhanced emission seen from 5 hr onwards. A similar time 
behaviour is visible both in the soft energy band (5=0.5- 



_ XSPEC model: wabs«mekal(Z=0.77 Z E ) 

EPIC (PN + M0S1+M0SZ) counts 
" S = 0.5-2.B keV 
- H=2. 8-8.0 keV 



log{N H /cm- 2 ) 




-0.5 0.5 

Hardness Ratio = (H-S)/(H+S) 

Fig. 3. Comparison between the observed hardness ratios 
and absorbed one temperature plasma models. The left 
and right hand vertical grey stripes indicate the hardness 
ratio values of the first and the second halves of the obser- 
vation, respectively, measured with EPIC (see lower panel 
of Fig. EJ. The curves give plasma temperature versus 
hardness ratio for a fixed value of the hydrogen column 
density. The dashed areas show the constraints on the 
model parameters when combining the results from the 
quantile analysis of PN data (see lower panel of Fig. 
with hardness ratios measured with EPIC. 



2.8 keV) and in the hard energy band (if=2.8-8keV). 
To quantify possible softening or hardening of the source 
events with time, we use the standard hardness ratio 
method by computing HR = (H — S) / (H + 5) (where the 
hardness ratio errors are computed using Gaussian error 
propagation) . The hardness ratios for the first and the sec- 
ond halves of the observation are HR = —0.29 ± 0.10 and 
0.12 ± 0.03, respectively (see the lower panel of FigureEJ- 
A least-square fit assuming constant hardness ratio 
leads to \ 2 — 15-4 (with one degree of freedom), with 
a probability 4 of Q ~ 10 -4 . This fit therefore excludes, at 
the 99.99% confidence level, the possibility that the hard- 
ness ratio remained constant during the XMM-Newton 
observation. We conclude that the X-ray emission from 
V16470ri hardened during the second half of the obser- 
vation. 

However during the second half of the observation, 
there is no correlation between the observed variations 
of the full energy band light curve and the hardness ratio. 
The hardness ratio values computed within time bins of 
1000 s (see the lower panel of Figure|2J) are consistent with 
a constant level (with a probability of 89%). In particular, 
the rise and the decay of the light curve do not correspond 
to any hardening and softening of the spectrum, whereas 
such a correspondence is usually observed in typical X- 
ray flares from young low-mass stars that are associated 



4 Here, Q is the probability that the best-fit model matches 
the data, given the value of % 2 ■ Values of Q < 0.05 indicate 
that the best-fit model is a poor fit, while Q ~ 1 indicates the 
fit is overconstrained. Optimal values of Q are ~ 0.5. 
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Fig. 4. Quantile analysis of the energy of the X-ray pho- 
tons of V16470ri collected by EPIC/PN. The upper panel 
shows the EPIC/PN light curve of V1647 0ri in the 0.5- 
8.0 keV energy band. Vertical dashed lines define variable- 
length time bins numbered from #1 to #5 where PN col- 
lected 200 net counts. For each time bin, the box plot 
indicates the energy distribution of the source (the ver- 
tical line and the box indicate the energy range and 
the quartiles/median, respectively) corrected from back- 
ground (except minimum and maximum values) using the 
Hong et al. (2004) method. The grey box plot correspond- 
ing to bin #2-#5 indicates the energy distribution of the 
entire second half of the observation. The lower panel 
shows the X-ray colour-colour diagram based on median 
and the ratio of two quartiles (Hong et al. 2004). The 
numbered asterisks mark the X-ray colours correspond- 
ing to the time bins defined in the upper panel. The di- 
amond marks the X-ray colours of the second half of the 
observation. Typical one sigma error bars are shown for 
comparison. 



with the heating and the cooling phases of the emitting 
plasma (e.g., ITsuboi et al. 19981 llmanishi et al. 2003f) : in 
such cases, the absorbing column density is observed to 
remain constant. 

Figure [3] shows a comparison of the hardness ratio val- 
ues measured with EPIC for the first and the second halves 
of the observation with absorbed one temperature plasma 



models computed with XSPEC (version 11. 3. 0) 5 . We can- 
not disentangle possible changes in absorbing column den- 
sity from changes in plasma temperature with only one 
hardness ratio. The usual improvement would be to intro- 
duce three consecutive energy bands, which would allow 
the computation of two hardness ratios: the first built on 
soft and medium energy bands, sensitive to the absorbing 
column density; the second built on medium and hard en- 
ergy bands, sensitive to the plasma temperature. However 
the total number of counts is too limited to permit us to 
obtain two hardness ratios with small enough error bars 
to allow a temporal study. To improve our analysis, we 
investigate changes in the energy of the source X-ray pho- 
tons collected by PN from quantile analysis, as proposed 



by Hong et al. (2004) 



To avoid signal-to-noise variations, we define variable- 
length time bins in which the total number of net counts 
collected by PN is constant, and equal to 200 net counts 
per time bin. This leads to 5 time bins(where bin #1 corre- 
sponds to the first half of the obser vation) and thereby al- 
lows a variability study. Following Hong et al. (2004)] we 
determine for each time bin median (E 50 %) and quartile 
(i?25%, £75%) energies of the source+background X-ray 
photons collected by PN, corrected for the background 
X-ray photons (see upper panel of Fig. ^J; and associ- 
ated quantiles, Q x % = (E x % - Ei ow )/(E up - E iow ), with 
Ei ow = 0.5 keV and E up = 8keV. Then, we compute 



the two X-ray colours proposed by Hong et al. (2004) 
x ee log(<2 5 o%/(l-Q 5 o%)), and y = 3 x Q 25% /Q 75% . These 
definitions of X-ray colours lead to a unique correspon- 
dance between the X-ray colour pair (x,y) and plasma 
parameters (N^,T), where iVu is the absorbing hydro- 
gen column density, and T the plasma temperature. The 
lower panel of Fig. 0] shows the reference mapping of the 
x,y space-parameter that we obtained by simulating with 
XSPEC absorbed one-temperature plasma spectra from 
an on-axis source observed with PN. This grid allows a 
direct estimate of logTVn and kT from x and y. Based 
on 10,000 trials of synthetic spectra computed from a 
model with parameters logTVn = 22.5, kT — 3.0 keV, 
Z = 0.77 Z Q (see g|and model #1 in Tabled), and 200 
counts, we found that the errors on x and y colours are 
distributed normally. In the lower panel of Fig. 01 we plot 
one sigma error bars on x and y for a single (200 net count) 
time bin. 

For bin taking into account the one sigma un- 
certainties on the X-ray colours, we find logiVn ~ 22.4 
(22.25-22.5), and kT ~ 2.2keV (1.8-2.8 keV). Over the 
course of the entire observation, the y colour varied by 
less than the one sigma uncertainty for a single time bin. 
From bin #1 to bins #3 and #4, there is a variation of the 
x colour greater than the one sigma error bar. However, 
due to the large error on the y colour, this variation in 



5 XSPEC wabsxmekal model, with global metallicity abun- 
dance Z — 0.77 Zq (obtained from spectral fitting; using 
the PN, MOS1, and MOS2 redistribution matrix files and ancil- 
lary response files appropriate for V1647 Ori (on-axis source) . 
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x colour implies no significant variation of the hydrogen 
column density and/or the plasma temperature. During 
the second half of the observation, we find logiVn ~ 22.53 
(22.49-22.57) and kT ~ 3.2 keV (2.9-3.7 keV); because 
~ 80% of the photons were detected in this latter half of 
the observation, the uncertainties on the quantile values 
obtained for the second half of the observation are smaller 
by a factor ~ 2 than the uncertainties in the first half of 
the observation. Therefore, there appears to be a hint of 
an increase in plasma temperature during the second half 
of the observation. 

In Figure |3 we combine these constraints on the hy- 
drogen column density and the plasma temperature of 
the first and the second halves of the observation with 
the hardness ratio values obtained previously. We con- 
clude that the hardening of the spectrum can be explained 
by a rather small increase of the plasma temperature 
from ~2keV to ~3keV, combined with a small increase 
of the absorbing column density from logTVn ~ 22.3 to 
logiVn ~ 22.5. 

It appears overall that the large increase of the count 
rate is not associated with any large variations of the tem- 
perature. This lack of strong correlation between X-ray 
flux and temperature suggests that we are perhaps also 
observing variations in the emission measure of the X- 
ray-emitting plasma. The variability observed with XMM- 
Newton and that measured previously with Chandra on 
2004 March 7 and 22 will be compared further in §5\ to 
investigate this possibility. 

4. Spectral analysis 

We created source and background spectra with asso- 
ciated redistribution matrix files and ancillary response 
files using the SAS task especget applied on the same 
source and background extraction region used for the light 
curve, and we binned the spectra to 15 counts per spec- 
tral bin using the ftools task grppha. Fig. [5] shows PN, 
MOS1, and MOS2 spectra of the whole observation af- 
ter suppression of high flaring background periods. The 
spectrum is clearly absorbed, with only a few counts de- 
tected by PN below 1 keV in energy, whereas PN detected 
counts up to 9 keV in energy. A prominent line is visible in 
PN spectra around 6.7 keV, corresponding to the FeXXV 
triplet emission line. We performed spectral modeling 
with XSPEC (version 11.3.0) using a multi-temperature 
plasma (mekal) model (Kaas tra et al. 1996JI with inter- 
vening absorption (as represented by the wabs model, 
which uses IMorrison fc McCammon 19831 cross-sections) . 
Table^gives the best fit parameters that we found for each 
statistically acceptable model (i.e. each model leading to a 
fit with Q- value greater than 5%), with the corresponding 
90% confidence level intervals. 

The absorbed one-temperature plasma model 
(wabsxmekal) leads to N H = 2.9 x 10 22 citT 2 (2.5- 
3.4xl0 22 cm-2), kT = 3.0 keV (2.4-3.8keV), with 
Z = 0.8 Z© (0.4-1.2 Zq), i.e. column density and tem- 
perature values consistent with the ones previously 
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Fig. 5. XMM-Newton X-ray spectra and our best fit mod- 
els of V1647 Ori. The upper panel shows the data binned 
to 15 counts in each spectral bin, and our best fit models 
with an absorbed two-temperature optically thin plasma 
(see model #2 in Table Continuous line, dashed line 
plus squares, and dotted line plus disks, stand for PN, 
MOS1, and MOS2, respectively. The lower panel shows 
the residuals of the fit, in units of the uncertainties in 
the individual data points. The FeXXV emission line at 
6.7 keV is clearly visible. 
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Fig. 6. Fit around the 6.7keV iron line region. The left 
hand panel shows an enlargement of the absorbed two- 
temperature model (see model 2 in Table around the 
6.7keV iron line, using a linear energy scale. The right 
hand panel shows the same energy area adding a Gaussian 
line fitting the 6.4 keV neutral iron line (see model 3 in 
Table Symbols are identical to the ones of Fig. [5] 

determined from our quantile analysis in J|l This quantile 
analysis remains the only time-dependent spectral analy- 
sis possible with this data set due to the small number of 
total counts. 

We obtain a fit with a better probability (Q=56%) 
when we add a second temperature component. Therefore, 
we consider this model (#2 in Table^) as our best fit. The 
bulk (~ 75%) of the emission measure in the 0.5-8 keV en- 
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Table 1. Best fit parameters of EPfC spectrum models with their errors at the 90% confidence levels (A% 2 = 2.71; 
corresponding to a = 1.64 for Gaussian statistics). For PN, MOS1, MOS2, the effective exposures after suppression of 
time intervals with high levels of flaring background are 33.2, 37.7, and 38.1 ks, and there are in the source extraction 
region 1264, 408, 408 counts, respectively (~ 20% of these counts come from the background). 



Temperature Emission Measure a 6.4 keV Fe Line b Lx,intr. a 

Model JVh.22 Z T soft T hard EM^ lt EM^ d N ph ,. 7 EW X l W° Q d 0.5-2.8/2.8-8/0.5-8 

# [cm" 2 ] [Z Q ] [keV] [cm" 3 ] [cm^s^ 1 ] [eV] [%] log [erg s" 1 ] 

1 2.9±g;| 0.8±g | 3-0lg-j 0.7±0.1 1.11 (129) 19 30.8 30.5 31.0 

2 4.1±0.6 0.8io;3 0.9±0.2 4.3±™ 1-2±J;1 0A±° o i 0.97 (127) 56 31.4 30.5 31.5 

3 4.1±0.6 0.8^3 0.9±0.2 i.2±*i l.Stli 0.4t£? i-^ttl 109 0.96 (126) 60 31.4 30.5 31.5 

a EM and Lx.intr. were computed assuming 400 pc for the distance of V1647 Ori. The X-ray intrinsic luminosity was estimated 

from the fit of the PN data. 
b Energy and width fixed to 6.4 keV and 12 eV, respectively. 

c v is the degree of freedom, and xt i s the reduced Chi-square (i.e. the Chi-square divided by the degree of freedom). 

The Q-value is the probability that one would observe the Chi-square value, or a larger value, if the assumed model is true, 
and the best-fit model parameters are the true parameter values. 



ergy band, and therefore most of the intrinsic X-ray emis- 
sion, arises with a soft plasma component, with kT so f t = 
0.9 keV (0.7-1.1 keV); whereas a hard plasma component, 
with fcThard = 4.2 keV (3.0-6.5 keV) contributes ~ 25% of 
the total X-ray emission. For this two-temperature model 
we find a column density a bit larger than in model #1, 
with N K = 4.1 x 10 22 cm~ 2 ( 3.5-4.7 x 10 22 cm- 2 ). 

One or two-temperature plasma models with the hy- 
drogen column density fixed to the lowest value found 
from optical/IR observations of V1647 Ori during outburst 
(e.g., Ay =11 mag from IVacca et al. 20041 see discussion 
below) lead to poor fits (Q < 0.05) and, hence, are not 
acceptable. 

The Kcv fluorescent line from neutral to low-ionized 
iron at 6.4 keV was detected by |Imanishi et al. (2001)| with 
Chandra from a class I protostar in the p Ophiuchi dark 
cloud during a bright flare, and was interpreted as re- 
flection emission from a face-on accretion disk. We esti- 
mate the 90% confidence limit on the strength of a neutral 
6.4 keV iron line in the V1647 0ri spectrum by adding a 
Gaussian line of 12 eV width at 6.4 keV. This will be im- 
portant to constrain future 2D models involving the heat- 
ing of plasma and its interaction with nearby circumstel- 
lar material. We retrieve exactly the best-fit parameters 
found for model #2 (see Table QJ, and the Gaussian line 
at 6.4 keV appears to improve the fit in the 6.7 keV iron 
line region (see Fig. EJ). However, we find that the total 
photon rate in the line is iV p h = 4.3 x 10 cm s _1 (0- 
8.7x 10 -7 cm -2 s _1 ) so, given this range of uncertainty, we 
cannot exclude the possibility that there is no 6.4 keV iron 
line in the XMM-Newton spectrum of V1647 0ri . 

Both the temperature and the column density found 
here with the one-temperature plasma model are con- 
sistent with results obtained from the same spectral 
modeling of Chandra data obtained on 2004 March 7. 
The large contribution from soft plasma emission to 
the X-ray spectrum of V1647 0ri revealed by XMM- 



Newton data is reminiscent of that observed from TW Hya 
l|Kastner et al. 2002jl . For TWHya, the soft X-ray emis- 
sion is believed to be generated via mass accretion from its 
circumstellar disk, which produces shocks onto the stel- 
lar photosphere. However the hard plasma emission of 
V16470ri above IkeV is too elevated to be explained 
in terms of accretion shocks onto a low-mass star, and 
can be understood only in the framework of plasma 
heating sustained by magnetic reconnection events (e.g., 
IKastner et al. 2004jl . 

This XMM-Newton measurement of the hydrogen col- 
umn density, consistent with the previous Chandra es- 
timate, is now well enough constrained to reveal a 
significant excess of hydrogen column density com- 
pared to the value derived from optical/IR observa- 
tions. Combining Nn/Aj = 5.6 ± 0.4 x 10 21 cm~ 2 mag" 1 
UVuong et al. 2003), and Aj/A v = 0.4008 - 0.1187 x 
(iJv/3.1)" 1 (|Cardelli et al. 1989|l . leads to Ay = 
(7V H /10 21 )/[2.24 - 0.66 x (Ev/3.1)" 1 ], where the param- 
eter Ry ranges from the average galactic value of 3.1 
to ^5.5 for a dusty environment. Therefore our mea- 
surement, iVn = 3.5-4.7 x 10 22 cm~ 2 , is equivalent to 
Ay =22.3-29. 6mag and 18.8-25.0mag, for R Y =3.1 and 
5.5, respectively. The highest value found for the visual 
extinction of V1647 0ri in its outburst is Ay — 11 mag 
i|Vacca et al. 2004jl . which is well below our 90% confi- 
dence interval. We stress that there remains a discrep- 
ancy between our lower range estimate and this high- 
est visual extinction for any value of the parameter Ry. 
We conclude that there is an excess of hydrogen col- 
umn density in the local environment of V1647 Ori of at 
least 1.5-2.6x 10 22 cm -2 , assuming the visual extinction 
of Vacca et al. (2004")] and Ry ~ 5.5. The gas-to-dust ra- 
tio along the line-of-sight towards V1647 0ri is hence un- 
usual, with more gas than dust towards this object. This 
is consistent with the picture of the rise of a wind/jet 
from the central embedded source, which would increase 
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Fig. 7. Long term multiwavelength variation of the observed flux density of V1647 Ori. The left hand panel (adapted 
from Kastner et al. 2004) shows the variation of the observed flux density of V1647 0ri from the optical (/ band), 
near-IR (J, H, and K bands), and X-ray (0.5-8 keV energy band). The first three X-ray measurements are from 
Chandra (Kastner et al. 2004), where we have computed for the 2002 November 14 measurement the 90% confidence 
level from observed counts using Gehrels (1986) statistics. The dotted arrow marks the X-ray upper limit at the 90% 
confidence level computed from the XMM-Newton observation of HH24-26 on 2003 September 3 (Ozawa et al., in 
preparation). The range of X-ray flux density that we report here with XMM-Newton is represented by a segment, 
detailed in the right hand panel. The right hand panel shows the X-ray flux density observed with XMM-Newton, 
computed from the EPIC light curve (upper panel of Fig. 0) using a conversion factor from count rate to flux density 
derived from our best spectral fit (Fig.[3J). One sigma error bars are given for both XMM-Newton and Chandra fluxes. 



the amount of gas close to central object where the X-ray 
emission is believed to arise. 

5. Comparison of V1647 Ori X-ray fluxes observed 
with Chandra and XMM-Newton 



We have updated Fig. 1 of Kastner et al. (2004) by adding 
IR measurements of Walter et al. (2004) which completes 
the monitoring of V1647 0ri through 2004 April, and 
XMM-Newton measurements (see the left hand panel of 
Fig. 0). XMM-Newton observed serendipitously the re- 
gion of V1647 0ri on 2003 September 3, with ~50ks expo- 
sure, by pointing toward the Herbig-Haro emission com- 
plex HH24-26 located about 4' south (Ozawa et al., in 
preparation). V1647 Ori was not detected. We convert this 
non-detection (computed on a ~25 ks exposure without 
flaring background) to a flux density upper limit assum- 
ing an absorbed (Nh = 5.7 x 10 22 cm~ 2 ) one-temperature 
(fcT=4.8keV) model with Z=0.3Z Q , i.e. the same model 
used by Kastner et al. (2004) to derive flux density from 
the Chandra detection of V1647 Ori on 2002 November 14. 
The resulting XMM-Newton upper limit is consistent with 
the low pre-outburst activity level of V16470ri observed 
by Chandra on 2002 November 14 during a 55.9 ks expo- 
sure. This demonstrates that the rise of X-ray activity had 



not yet occured two months before the optical/IR out- 
burst and strengthens the apparent correlation between 
the sharp increase in X-ray flux from pre-outburst to out- 
burst and the similar brightening seen in the optical/IR, 
already noticed by Kastner et al. (2004) Using our best- 



fit spectral model, we convert EPIC count rate in the 0.5- 
8keV energy range (see Fig. P3J) to observed X-ray flux 
density (note that the conversion factor is essentially inde- 
pendent of the number of plasma components). As shown 
in the right hand panel of Fig. the X-ray flux observed 
with XMM-Newton ranges from roughly the flux observed 
by Chandra on 2004 March 22 (i.e. ~10 times greater than 
the pre-outburst X-ray flux) to a peak at a value two times 
greater than the one caught by Chandra on 2004 March 7 
(^200 times greater than the pre-outburst X-ray flux). 
The X-ray variability of V1647 0ri in outburst is hence 



also clearly enhanced , as noted by Kastner et al. (2004) 



The right hand panel of Fig. reveals that the X- 
ray flux density of V1647 0ri during the first half of the 
XMM-Newton observing period on 2004 April 4 is nearly 
identical to the flux density of the source during the epoch 
of the Chandra observations on 2004 March 22. Similarly, 
the enhanced flux density of V1647 0ri during the second 
half of the XMM-Newton observing period is an excel- 



Grosso et al.: Enhanced X-ray variability from V1647 0ri 



Willi I' ''Ml i 



. ; ■: i ■ . 

|l|l IjjlJ, | I I 



Trial period [day] 




Fig. 8. Folded XMM-Newton and Chandra X-ray fluxes of V1647 Ori. The left hand panel plots the function V 
versus trial period, from the value given by the length of the XMM-Newton observation to 2.3 days with period step of 
10 -6 s. This function is defined only for periods allowing the folding of the two Chandra measurements with the XMM- 
Newton light curve (see text for details) . Diamonds mark the best three shortest trial periods with their corresponding 
n th -harmonic. The right hand panel displays the resulting XMM-Newton (line) and Chandra (grey hourglasses) flux 
measurements folded using these three periods. The error bars represent the equivalent XMM-Newton flux integrated 
during the Chandra observing time. The last plot of the right hand panel shows for comparison a bad fit corresponding 
to the point marked by an asterisk in the left hand panel. 



lent match to the X-ray brightness of the source during 
the epoch of the Chandra observations on 2004 March 7. 
The combined Chandra+ XMM-Newton data set strongly 
suggests the possibility that the V1647 Ori shows a peri- 
odic variation in X-ray brightness. We have investigated 
this possibility by folding the Chandra flux measurements 
with those from XMM-Newton. 

For this purpose, we define the following function: 



V(P) 



\ 



y {Fcxo 

i=l 



FxMM,i,p) 2 



a XMM,i,P 



where P is the trial period; Fcxo,i is one of the two 
Chandra flux measurements; Fxmm,%,p is the XMM- 
Newton flux integrated over the Chandra observational 
window (depending of the trial period) using a linear in- 
terpolation of the XMM-Newton data points; <JxMM.i,p 
is the standard deviation of the XMM-Newton flux in 
the Chandra observational window. This function is de- 
fined only for periods allowing the folding of the two 
Chandra measurements with the XMM-Newton light 
curve. T> defines in the folded light curve the quadrat- 
ically averaged distance between Chandra measurements 
and the XMM-Newton light curve, in units of the standard 
deviation of the XMM-Newton flux in the Chandra obser- 



vational window. We computed numerically T> versus trial 
periods, from the value given by the length of the XMM- 
Newton observation to 2.3 days, with period step of 10~ 6 s. 
The left hand panel of Fig. |H| shows the resulting value of 
T> for each trial period. The structure of the T> function 
thus obtained is non continuous and cyclical, thereby ap- 
pearing similar to an attractor. We consider the periods 
providing the smallest value of T> as the best period candi- 
dates. The first three period candidates are: P=0.43, 0.54, 
and 0.72 day, with V ~ 0.008. The other periods leading 
to a minimum value for T) are n th -harmonics (where n is 
an integer) of these first periods. The right hand panel of 
Fig. |21 shows the corresponding X-ray folded light curves 
obtained with these period candidates. The last plot in 
Fig. El illustrates a typical period- folding result when V is 
greater than 1, corresponding to the asterisk point in the 
left hand panel of Fig. |SJ clearly, in this case, we obtain 
a poor fit. Because our combined Chandra and XMM- 
Newton dataset spans nearly 1 month, any folding of the 
Chandra flux measurements with periods much smaller 
than ~1 month is very sensitive to the assumed value of 
the period. 

From both an observational and a physical point 
of view, the 0.72 day period appears the most likely. 
Observationally, with this period, roughly 60% of the 
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phase is covered by the XMM-Newton observation, where 
during half of this time interval the X-ray flux is close to 
one of the two levels measured with Chandra. Of course if 
we assume that the light curve is indeed symmetric about 
phase 0.5, then the X-ray flux is, during the whole period, 
close to one of the two levels measured with Chandra pre- 
cisely half the time. Therefore random short X-ray obser- 
vations would catch easily one of these two characteristic 
level of activity, explaining naturally the results obtained 
with Chandra. Physically, this period corresponds to the 
time scale of Kcplerian rotation at a distance of 3.4 so- 
lar radius from a one solar mass star, or equivalently to a 
distance of 1 (1.4) stellar radius for a one solar mass star 
with an age of 0.5 (1.0) Myr IjPalla fc Stahler 1999jl . 

6. Discussion 

Until very recently, no X-ray satellites have had the oppor- 
tunity to observe a PMS eruptive object at the beginning 
of its outburst. In this respect the recent Chandra and 
XMM-Newton observations of V16470ri offer new insight 
into a PMS accretion burst and, hence, into PMS ac- 
cretion processes. During phases of normal accretion ac- 
tivity in T Tauri stars (Af acc ~ 10" 7 M Q yr~ 1 ) the dy- 
namic pressure from the disk accretion flow is balanced 
by the magnetic pressure of the stellar magnetosphere at 
the Alfven radius l|Ghosh fc Lamb 19781 |Koenigl 19911 ). 
The magnetic pressure inside this radius (close to the co- 
rotation radius) is strong enough to prevent the accreting 
material in the disk midplane from spiralling towards the 
star, but rather lifts it into free fall along the stellar field 
lines, leading to magnetospheric accretion (see review by 
IShu et al. 20001) . 

In contrast, during an accretion outburst, the ram 
pressure from the accretion disk flow is sufficient to 
push the inner edge of the accretion disk within the co- 
rotation radius, leading to strong magnetic interaction 
between the star magnetosphere and its inner accretion 
disk. This interaction area, called the reconnection ring, 
is then a possible source of hard and soft X-ray emis- 
sion from magnetic reconnection events (|Shu et al. 19 97). 
Strong reconnection events will also occur when differ- 
ential rotation between the inner disk edge and stel- 
lar magnetospshere cause fields lines to twist, strengthen 
and expand l|Lovelace et al. 19951 IGoodson et al. 19971 
IMatt et al. 2 002). In these cases reconnection drives out- 
flows occuring periodically above and below the disk 
( IGoodson fc Winglee 1999| IGoodson et al. 1999|l . 

In the most extreme case, the accretion flow is even suf- 
ficient to crush the stellar magnetic field against the stellar 
photosphere, leading to the formation of a boundary layer 
at the inner edge of the accretion disk ( Hartmann 1998), 
which may quench the star-disk X-ray emission. 

The previous shorter (5 ks exposures) observations of 
V1647 0ri obtained with Chandra have shown on 2004 
March 7 during the outburst an enhanced X-ray flux com- 
pared to the pre-outburst level measured by Chandra on 
2002 November 14; and a drop of the X-ray flux on 2004 



March 22, which we previously interpreted as the possible 
onset of a quenching X-ray emission phase, or the trigger- 
ing of a phase of strong variability in both X-ray lumi- 
nosity and temperature (JKastner et al. 2004). The longer 
~ 38 ks observation with XMM-Newton of V1647 0ri on 
2004 April 7, presented here, shows that there is indeed 
as yet no quenching X-ray emission phase, but there is 
clearly enhanced X-ray variability. 

The detailed X-ray spectrum obtained with XMM- 
Newton improves of our knowledge of the X-ray properties 
during the outburst phase, showing that ~ 75% of the in- 
trinsic X-ray emission in the 0.5-8 keV energy band may 
be explained by an accretion shock onto the photosphere 
of this low-mass star, whereas ~ 25% of the total X-ray 
emission comes from a plasma heated by a magnetic re- 
connection mecanism. The dramatic increase of accretion 
rate, producing the optical/IR outburst, may also explain 
the observed rise of X-ray emission from V1647 Ori: first 
by increasing the quantity of gas falling onto the stellar 
surface, which radiates soft X-rays at the accretion shock; 
second by pushing the inner accretion disk boundary in- 
side the corotation radius, which would increase the con- 
tribution of the reconnection ring to the soft/hard X-ray 
emission. 

Our quantilc analysis has shown that large increase 
of the count rate is not associated with any large tem- 
perature variations as for typical X-ray flares from YSOs, 
suggesting that we are likely also observing variations of 
the emission measure. The variation of the emission mea- 
sure can be produced by a variation of the (observed) 
plasma volume and/or a variation of the plasma electronic 
density. We cannot totally exclude that we observed an 
unusual X-ray flare with a strong variation with time of 
the plasma electronic density in the flaring loop. However, 
the comparison of our XMM-Newton light curve with the 
Chandra measurements has suggested that the X-ray flux 
of V1647 0ri could be periodic, with a likely period of 
0.72 day - corresponding to the Keplerian rotation of close 
material in the close vicinity of the star. Therefore, the 
enhanced X-ray variability of V1647 Ori is more likely ex- 
plained by a variation of the observed plasma volume by 
rotational modulation. In this scenario, the inner part 
of the accretion disk eclipses periodically a fraction of 
the plasma volume close to the stellar surface, produc- 
ing a dimming of the observed X-ray flux by a factor 
of about 20. Rotational modulation by accretion funnels 
has been proposed to explain optical eclipses observed 
in the classical T Tauri star AATau, where the accre- 
tion disk is thought to be warped by the magnetospheric 
magnetic dipole inclined respectively to the star's spin 
IJBouvier et al. 199911. 

We stress that the optical/IR flux of V1647 Ori is dom- 
inated by the outburst luminosity of its inner accretion 
disk and scattered light from the surrounding nebulosity, 
therefore any optical/IR modulations of the photospheric 
flux are difficult if not impossible to observe. Therefore X- 
rays, likely produced only close to the central star, offer a 
unique observational window to probe the star-disk inter- 
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action area during this PMS evolutionary phase in which 
the accretion rate, and hence the source luminosity, are 
dramatically elevated above normal. 

Longer X-ray observations with XMM-Newton and 
Chandra will be needed to confirm both the candidate pe- 
riod, and to validate the scenario proposed here to explain 
the enhanced X-ray variability. 
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